Elastic scattering of a quantum matter-wave bright soliton on a barrier 
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We consider a one-dimensional matter-wave bright soliton, corresponding to the ground bound 
state of N particles of mass m having a binary attractive delta potential interaction on the open 
line. For a full iV-body quantum treatment, we derive several results for the scattering of this 
quantum soliton on a short-range, bounded from below, external potential, restricting to the low 
energy, elastic regime where the centre-of-mass kinetic energy of the incoming soliton is lower than 
the internal energy gap of the soliton, that is the minimal energy required to extract particles from 
the soliton. 

PACS numbers: 03.75.Gg, 03.75.Lm, 34.50.-s 



I. INTRODUCTION 

The ultracold atomic Bose gases have proved to be very 
flexible physical systems, where both the dimensionality 
and the interaction strength can be adjusted at will. By 
trapping Bose-condensed atoms in an optical waveguide 
that freezes their transverse motion in its ground state, 
one obtains ultracold one-dimensional Bose gases. By 
further making the effective one-dimensional atomic in- 
teraction attractive, one can produce matter- wave bright 
solitons, which are bound states of matter with typically 
thousands of particles (H-Q- This opens up a new field 
with exciting possibilities, the field of coherent matter- 
wave optics with massive objects. Even a soliton of light 
atoms such as 7 Li is typically more massive than the big 
organic molecules (such as fullerenes) used in interfero- 
metric experiments *4] , and it has the advantage of having 
much larger centre-of-mass de Broglie wavelength, since 
the atomic gases can be prepared in the nK temperature 
range 

The disadvantage (or depending on the perspective, 
the additional feature) of the matter- wave soliton is that 
it is a quite fragile object: The ground state soliton is sep- 
arated by a continuum of fragmented solitons by a small 
energy gap A, with A/fee typically sub-microKelvin. For 
the scattering of a quantum soliton on a barrier to be 
guaranteed to be elastic by energy conservation, one has 
to restrict the kinetic energy of the centre of mass of the 
soliton to values below the gap A. This elastic scattering 
regime is quite intriguing and was recently considered in 
proposals for production of real space Schrodinger-cat- 
like states by coherent splitting by a laser barrier of the 
centre-of-mass wavepacket into transmitted and reflected 
components @, 0] and for Anderson localisation of quan- 
tum solitons in a disordered potential [§[ [33j]. On the 
experimental side, scattering of a soliton on a barrier is 
under experimental investigation, for the moment out of 
the elastic scattering regime, with fragmentation of the 



soliton into two main pieces [Tol | . 

Here we restrict to the elastic scattering regime on a lo- 
calized potential barrier: inside or close to the potential, 
the system can virtually access internal excited states 
(where the soliton is fragmented) but it fully occupies 
the ground state soliton at asymptotically large distances 
from the barrier, so that scattering of the soliton with in- 
coming centre-of-mass wavevector K is characterized by 
the transmission amplitude t and the reflection ampli- 
tude r with |r| 2 + |i| 2 = l. As an initial wavepacket may 
be expanded over such stationary scattering states, its 
time-dependent wavefunction away from the barrier can 
be deduced from the A"-dependent t and r amplitudes. 

Whereas the classical field (or Gross-Pitaevskii) equa- 
tion was extensively used to study soliton dynamics and 
fragmentation in external potentials [111 ], it does not 
look appropriate in the elastic scattering regime. First, 
the Gross-Pitaevskii equation does not provide a full 
quantum-mechanical treatment of the centre-of-mass mo- 
tion. In the absence of an external potential, it predicts 
the existence of localized stationary solutions, whereas 
the centre-of-mass position necessarily spreads ballisti- 
cally in time in the quantum world |12| . In the scatter- 
ing by a barrier, it cannot describe Schrodinger-cat-like 
states, where the unfragmented soliton has some non- 
zero probability amplitude to be to the left (resp. to the 
right) of the barrier Q. Secondly, the classical soliton 
misses the rigidity of the quantum soliton at the heart of 
elastic scattering: in the classical field theory, the mov- 
ing soliton can in principle always slow down by radiating 
at infinity an arbitrary small amount of energy, without 
violating energy conservation, whereas in the quantum 
theory, the number of particles radiated to infinity (that 
carry away an energy at least A) is quantized. 

We thus have to use the quantum field theory, which 
constitutes a full many-body problem when the number 
of bosons N is large. In the absence of a barrier, it 
was solved with the Bethe ansatz generalized to com- 
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plex quasi- moment a, both for the ground state [13j and 
for the excited states fT^ - fl6| : the many- particle ground 
state in the presence of a harmonic trap was investigated 
in [T3|. In the presence of a barrier, the Bethe ansatz is 
not applicable, the exact iV-body solution is not known 
and one has to resort to approximations. When the bar- 
rier is broad as compared to the soliton size, it is natural 
to introduce the average V(X) of the external poten- 
tial experienced by the N bosons over the density profile 
p(x\X) of the ground state soliton with centre of mass 
localized in X. Then one writes a Schrodinger equation 
for a centre-of-mass wavefunction &(X), treated as a sin- 
gle particle of mass M — Nm (m is the mass of a single 
boson) moving in the potential V{X). This intuitive ap- 
proximation was used e.g. in [fl HJ. 

The scope of the present paper is to provide tools to 
construct this approximation, to control it with rigor- 
ous error bounds on the transmission and reflection am- 
plitudes, and to go one step beyond it in the large-TV 
limit. In Sec. [HI we define the problem; using a projector 
technique, we show that the centre-of-mass wavefunction 
can be given a precise meaning and that it obeys, 
in the elastic regime, an exact Schrbdinger-like equation 
with an effective potential that, in addition to V(X), 
contains a non-local and energy dependent contribution 
8V originating from all possible virtual fragmentations of 
the soliton. In Sec lIIII we derive a simple upper bound on 
the matrix elements of 8V, which allows to derive upper 
bounds (already used in |6j) on the error on t and r due 
to the omission of 8V; in the case of a very narrow poten- 
tial barrier, such as a repulsive Dirac delta, we show how 
to improve the procedure to get usable upper bounds. In 
Sec lIVI we determine from Bogoliubov theory the lead- 
ing order contribution to 8V in the large N limit, with 
this limit constructed in such a way that V(X) remains 
fixed. In Sec. [V] we again consider the large- N limit case 
with a Born-Oppenheimer-like approach, the heavy par- 
ticle being the centre of mass, and we identify a regime 
where it approximately coincides with the Bogoliubov re- 
sult of Sec lIVI In Sec lVIl we give simple applications of 
the formalism. We conclude in Sec. IVIII 



In first quantized form: 



H = H + V 
Ha 



N 



i=l l<i<j<N 
N 



V 



(i) 

(2) 
(3) 



where Xi is the spatial coordinate of the i th boson and Pi 
is its momentum operator. 

The free space Hamiltonian Hq can be diagonalised 
with the Bethe ansatz [T34l6j . Another key feature of 
Hq, that we shall use extensively, is the separability of 
the centre-of-mass degrees of freedom (associated to the 
centre-of-mass position X) from the internal degrees of 
freedom (whose N — 1 spatial coordinates can be ex- 
pressed in terms of the (xi)i<i<N through Jacobi for- 
mulas [19( , that are not required here) . This gives a ten- 
sorial product structure to the Hilbert space between the 
centre-of-mass variable and the internal variables, and it 
corresponds to the following splitting between the centre- 
of-mass kinetic energy operator and the internal Hamil- 
tonian H mt : 



Hn = 



P 2 

2M 



Hi, 



(4) 



with M = Nm is the total mass and P — Y\ pi is the to- 
tal momentum operator. The internal Hamiltonian H mt 
does not depend at all on the centre-of-mass variable [13] ■ 
It has only one discrete eigenstate, its ground state \<j>) 
of eigenenergy E (N) given by [l3[ 



1)JV(JV + 1). 



(5) 



Considered as a function of the Xi's, <f> is also the ground 
state of Hq since it corresponds to a centre of mass at 
rest, H \tf>) = Eq(N)\4>). It has a simple expression in 
terms of the single particle coordinates [l3[ 



, xn) = Me ^ ^i<i<]<_ 



(6) 



II. DEFINITION OF THE PROBLEM AND THE 
EFFECTIVE POTENTIAL FOR ELASTIC 
SCATTERING 

A. Hamiltonian and free space properties 

We consider N spinless bosons of mass m moving 
quantum-mechanically on the open one-dimensional line. 
The bosons have an attractive Dirac pair interaction 
characterized by the negative coupling constant g [l8j . 
and each boson is subjected to a localized potential U(x), 
that is U(x) rapidly tends to zero for |x| — > +oo. The N- 
body Hamiltonian H is the sum of the free space Hamil- 
tonian Hq and of the external potential Hamiltonian V. 



with a normalisation condition also easily expressed, by 
fixing the centre-of-mass position to the origin of coordi- 
nates: 



dxi . . 
which leads to [14|, 21 
\M\ 2 



dx N 8 ^y)xj /N^J \4>{xi,. . . , 



xn)\ 2 = 1, (7) 



m\g\ 



N-l 



N 



(N-iy. 



= i. 



(8) 



Apart from this discrete eigenstate, the spectrum of H- mt 
is a continuum separated from Eq(N) by a gap A, cor- 
responding to any possible fragmentation of the ground 
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state soliton into smaller solitons (including single parti- 
cles) with arbitrary centre-of-mass momenta. From the 
full spectrum obtained by the Bethe ansatz [ill - fill ], and 
using E (ni) + E {n 2 ) > E Q (ni + n 2 ) for n\,n% > 0, one 
finds 

A = E (N - 1) + E (l) - Eo(N) = ^N(N - 1), (9) 

i.e. A is the energy required to extract a particle of van- 
ishing relative momentum from the ./V-particle soliton. 

In the presence of the external localized potential U(x), 
the centre of mass of the gas experiences scattering and 
is no longer decoupled. Let us assume first that U (x) is 
everywhere non-negative, so that no boson can remain 
trapped in the potential. Then at low enough energy E 
such that 



E (N) <E< Eq(N) + A 
the eigenstate \ip) of H of energy E, 
= (E-H)\4>), 



(10) 

(11) 



has a simple structure far away from the external po- 
tential, corresponding to elastic scattering of the ground 
state soliton: Because of energy conservation, at a centre- 
of-mass position X — > ±oo, the internal state is in its 
ground state (f> and the centre-of-mass wavefunction as- 
sumes the usual asymptotic form of a single particle scat- 
tering state. Introducing the positive K such that 



, . h 2 K 2 , h 2 K 2 A 

E = E (N)-\ with < A. 

uv ' 2M 2M 

we thus have the boundary conditions (see Fig. [T|): 



(12) 



1p(x!,. . .,X N ) 

$P0 



X- 



$(X)<f>(x 1 ,...,x N ) (13) 
(14) 



JKX 



-iKX 



$(X) ~ te 

X^+oo 



iKX 



(15) 



where |r| 2 + |t| 2 = l. As we shall see, a meaning can 
be given to the centre-of-mass wavefunction &(X) of the 
soliton at all X, not simply at infinity. The goal of the 
present work is to calculate approximately the reflection 
amplitude r and the transmission amplitude t, and to 
control the resulting error. 

The previous physical reasoning has to be adapted 
when U (x) presents weakly negative parts, that may 
support bound states, in which case the scattering state 
could be fragmented, e.g. it could correspond to a (N — 
n)-particle soliton flying away, with n bosons trapped 
within the external potential (0 < n < N). A rt-particle 
bound state, being an eigenstate of the free space Hamil- 
tonian plus external potential Hamiltonian, has an en- 
ergy -Ebound necessarily larger than the sum of the mini- 
mal eigenvalues of each Hamiltonian, -©bound > Eq(ji) + 
ninfj, U(x). The energy of the fragmented scattering 
state is thus larger than Eo(N—n)+Eo(n)+ri inf x U(x) > 



Eq(N) + A + Nmi x U(x). Fragmented scattering states 
are thus forbidden by energy conservation over the energy 



range 



E (N) <E <E (N) + A + Nw£U(x), (16) 



a constraint over E and inf x U (x) that we assume to be 
satisfied (to guarantee purely elastic soliton scattering) 
and that will be recovered by a purely mathematical rea- 
soning. 



B. The effective potential 

An exact rewriting of Schrodinger's equation within a 
restricted subspace is obtained with the action of projec- 
tors on the resolvent G(z) = l/(z — H) of the Hamil- 
tonian [22j . leading to an effective Hamiltonian. Using 
the tensorial product structure of the Hilbert space be- 
tween the centre-of-mass and the internal variables, we 
define the operator projecting orthogonally the internal 
variables onto their ground state \<f>): 



V = 1cm® 10X01 



(17) 



where 1cm stands for the operator identity over the 
centre-of-mass variables. The supplementary orthogonal 
projector is 



Q = l-V. 



(18) 



Then for any complex and non-real number z, we obtain 
the exact expression [23j | 



VG(z)T 



V 



zV - VHV - VVQ 



Q 



zQ.-Q.HQ 



QVV 



(19) 



To access the scattering state of energy E, one should 
usually take the limit z — E + ie, e — > + , because the 
operator zQ — QHQ is usually not invertiblc for z = E 
(within the subspace over which Q projects). But wc 
shall now restrict to a situation where this operator is 
invertible because it is strictly negative. To this end, we 
assume that the external potential U(x) is bounded from 
below, 



inf U(x) > —00 



(20) 



Then the spectrum (abbreviated as Spec) of the Hermi- 
tian operator QHQ (within the subspace over which Q 
projects) is also bounded from below, 

inf Spec QHQ > N inf U(x) + E (N) + A. (21) 

X 

To ensure that QHQ — EQ is strictly positive (within 
the subspace over which Q projects), we thus impose 



h 2 K 2 

2M 



< A + TV inf U(x), 



(22) 
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which reproduces the physical result (ITS)) . The action of 
the projector V onto the eigenstate gives 



P\tf>) = |*) 



The wavefunction 



= (X\$) = ((X\ ® 



(23) 



(24) 



plays a crucial role, it is the centre-of-mass wavefunc- 
tion within the subspace where the internal variables are 
in their ground state (that is in the minimal energy N- 
particle soliton). In other words, &(X) is the centre-of- 
mass wavefunction of the soliton. In what follows, we 
shall use a shorthand notation of the type |$) = (<f>\ip), 
where the tensorial product structure between centre-of- 
mass and internal variables is implicitly assumed. In 
terms of the original variables xt , it is expressed as 

N 



$(Jf ) = J dx-L... dx N S(X - Xi/N) 



<fi*(xi,...,x N )il>(x 1 ,...,x N ). (25) 

For z — ¥ E, the effective Hamiltonian appearing in the 
denominator of Eq. (|19l) is Hermitian under condition 
([22]) . and we find the exact Schrodinger-like equation for 



h 2 K 2 

2M 



1$) = 



P 2 
2M 



+ V(X) + SV 



(26) 



where P = —ihdx and, as we shall discuss, V is given 
by (f29|) and SV is given by ff3T~|) . This result can also be 
obtained by a direct calculation without introducing the 
resolvent: One applies the projector V and the projector 
Q to Schrodinger's equation (fTT|) . and one inserts the 
closure relation V + Q = 1 to the right of (E — H) to 
obtain 



{EQ 



EP\t/>) = VHV\ip) 
QHQ)QW) = QHVty) 



(VHQ)Q\iP) (27) 
(28) 



Under the condition (|2"2"|) we multiply the second equa- 
tion by the inverse of the operator EQ — QHQ (within 
the subspace over which Q projects) and we report the 
resulting value of Q\ip) within the first equation. After 
some rewriting using in particular (|23j) we recover (|26p . 

The first contribution to the effective potential in (|2"6")l 
is very intuitive, and is a simple function of the centre- 
of-mass position X, 



V(X) 



\v\<t>) = (V) 



x 



(29) 



where we have introduced (. . .)x the expectation value 
in the internal ground state for a fixed value X of the 
centre-of-mass position. For a general observable O that 
is diagonal in terms of the original Xi variables, one has: 



N 

(0) x = J dxx... dx N 8(X - Xi/N) 



x 0(xi, . . .,x N )\(f>{xi,. . . ,x N )\ 2 . (30) 



We shall see that V(X) is then simply the convolution 
of the external potential U(x) with the density profile 
p(x\0) of the ground state soliton whose centroid is fixed 
at the origin of coordinates. 

The second contribution to the effective potential in 
Eq. (|26l) is both non-local and dependent on the scatter- 
ing state energy E: 



SV 



\VQ 



EQ-QHQ 



QV\4>). 



(31) 



Its evaluation cannot even by performed with the Bethe 
ansatz, due to the presence of the external potential V 
in the denominator. A first strategy is to simply neglect 
SV as compared to V in ([26| , as already done in @, 0] , 
which intuitively should be accurate in the large N limit 
and/or when U(x) is broad as compared to the soliton 
size £. In Sec. IIIIl rigorous bounds on the resulting er- 
ror on the soliton transmission and reflection amplitudes 
are given. A second strategy is to calculate the leading 
large- N asymptotic expression of SV within Bogoliubov 
theory, as done in Sec. IIV1 or to rely on a simpler Born- 
Oppenheimer-like approximation, as done in Sec.|Vl that 
can be subsequently used in a numerical solution of (|26l) . 



C. How to calculate V(X) 

To obtain an operational expression for V(X), we in- 
troduce the mean density of particles in the ground state 
soliton for a fixed position of the centre of mass, 

p(x\X) = (p(x)) x , (32) 

where the operator giving the density in point x is 

N 

p(x) = y^ 5(xj - x). (33) 

8=1 

Using Eq. (13T))) we thus reach the intuitive result 



V{X)= I dxU(x)p(x\X) = / dxU(X - x)p(x\0) 
Jn Jm 

(34) 

using the translational invariance p{x\X) = p(x — X\0) 
and the fact that p(x\0) is an even function of x. This 
also enables an exact evaluation of V(X), since p(x\X) 
was calculated with the Bethe ansatz in [3, [2l[ : 



N-2 



= jvf ^ (JV - 2 - fc)!(JV + fc)! ' 



(-l) fc (fc + l) B -(fc+i)l»- 



x\/i 



(35) 



where f is the spatial width of the classical field (Gross- 
Pitaevskii) soliton, 



m\g\N 



(36) 
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A large- N expansion can be obtained from (|35| (2ll . l2~ij : 
For AT — > oo with Afg (and thus £) fixed, 

p(x|0) = A^(x) - C a ^[^(»)] + (37) 

where the classical field soliton single particle wavefunc- 
tion (normalized to unity) is given by 



4> (x) 



1 



1 



(38) 



For N —, 
leads to 



2£V2 cosh(a;/2£)' 
-oo with fixed £, a double integration by part 



V(X) = / dxN(j)l{x) 



[/(.T + A) - ^[/"(x + X) 



o(U). (39) 



III. BRACKETING THE TRANSMISSION AND 
REFLECTION AMPLITUDES 

In the simplest treatment of elastic soliton scattering, 
one simply neglects the contribution SV in (f26| The 
challenge is to be able to put a bound on the correspond- 
ing error performed on the transmission and reflection co- 
efficients. To this end, we derive an upper bound on the 
modulus of the matrix elements of SV. Then we derive 
upper bounds on the contribution of SV to the transmis- 
sion and reflection coefficients, first in a minimal version 
(where the bounds can be directly evaluated from exist- 
ing Bethe ansatz results) , and second in a refined version 
applicable also to arbitrarily narrow external potentials 
(such as a repulsive delta potential). 



A. Upper bound on the matrix elements of SV 

Let us consider two kets |$i) and |$ 2 ) for the centre-of- 
mass degrees of freedom. The corresponding wavefunc- 
tions need not be square integrable but the following kets 
should be normalisable, 



|«i> = QV\$i) 
\ U2 ) = QV\$ 2 ) 



(40) 
(41) 



The matrix element of SV may thus be written as 

($!|5V|* 2 ) = (ui|(-0)|«2) (42) 
where we have introduced the Hermitian operator 

Q 



QHQ - EQ' 



(43) 



According to relations (|21I22[) . the operator Q is positive 
(within the subspace over which Q projects). From the 
Cauchy-Schwarz inequality, 



From Eq. (1211) we find, e.g. injecting a closure relation in 
the eigenbasis of QHQ, 



(u 2 \G 2 \u 2 ) < 



\u 2 \u 2 ) 



[Nin£ x U(x) + A - h 2 K 2 /2M} 2 
Calculating the norm squared of the we find 

where we have introduced the positive quantity 
w{X) = (V 2 ) x {V)\ 



(45) 



(46) 



(47) 



with (. . .)x defined in ([501) and V given by ([3]). 

In conclusion, the non-local contribution SV to the ef- 
fective Hamiltonian appearing in (|26l) can be bounded by 
the local positive potential (that we call error potential) 



W(X) 



w(X) 



N mi x U(x) + A - h 2 K 2 /2M 
in the following sense: 

|<*i|*v|$2>| < f($i|^(x)|$ 1 )($ 2 |vK(A > )|$ 2 ; 



(48) 



1/2 



(49) 



B. Upper bound on the error on the scattering 
coefficients 



In the simplest approximation, one neglects SV in 
Eq. (|26| and one calculates the wavefunction §o(X) cor- 
responding to scattering of the soliton centre of mass onto 
the potential V(X) with incoming wavevector K > (i.e. 
for a soliton coming from the left): 



h 2 K 2 

2M 



■MX) 



2M 



*%{X) + V{X)* (X). (50) 



It leads to a transmission amplitude t$ and a reflection 
amplitude r , see Fig. [1] 

In the exact treatment, keeping the transmission 
and reflection amplitudes are t and r, see Fig. [TJ We 
introduce the two positive quantities: 



V = 



M($ \W(X)\$o) 

h 2 K\t \ 

($ \W(X)\Q ) 
($ Q \W(X)\$ ) 



(51) 
(52) 



where the scattering solution onto V(X) for an incoming 
wave with negative wavevector —K (i.e. for a soliton com- 
ing from the right) can be expressed in terms of <5>q(X): 



MX) 



$5(A)-r *$ (X) 



(53) 



\(ui\g\u 2 )\ 2 < (ui| U i)( U2 |g 2 | U2 ). 



(44) Then we have the rigorous result: 



6 



V(X)+5V 



V(X) 



*(X): 



3> n (X) 




FIG. 1: Left panel: In the regime of elastic scattering of a 
iV-boson bright quantum soliton on a localized external po- 
tential, the centre-of-mass wavefunction &(X) obeys bound- 
ary conditions for X — > ±00 with transmission and reflection 
amplitudes t and r satisfying \t\ 2 + \r\ 2 = 1. The effective po- 
tential seen by $ involves a non-local part SV corresponding 
to virtual soliton fragmentation, see Eq. ((26} . Right panel: In 
the simplest approximation, one neglects the non-local part 
SV, only keeping the average V(X) of the single particle ex- 
ternal potential U(x) over the density profile p(x\X) of the 
soliton with centre-of-mass in X. The resulting approximate 
transmission and reflection amplitudes are to and ro. Bounds 
on |t — to I and \r — ro| are derived in Sec. Mil and tools to 
evaluate the leading order corrections due to SV in the large 
N limit are developed in Sec. IIVI and in Sec. fVl 



Theorem: If eq I 2 < 1/2 then 

||*| - |*o 1 1 < I* — * | < 
IM - Ml < \r - r \ < 



\to\er) 



1/2 



1 - 2er? 1 /2 

l*ok 
1 - 2e?7 1 /2- 



(54) 
(55) 



The first inequalities result from the triangular inequal- 
ity. The proof of the second inequalities is given in the 
Appendix [A] In the case of an even external potential, 
one has simplified relations, since r\ = 1: 
Theorem: IfU(x) is even and e < 1/2, then 



||*|-|to||<|*-*o| < 
\r\ - \r \\ <\r- r \ < 



1 - 2e 
\to\e 
1 -2e' 



(56) 
(57) 



This writing immediately reveals that w(X) depends on 
correlations that would not be accurately treated in the 
classical field (Gross-Pitaevskii) approximation. From 
the Bethe ansatz wavefunction © of the soliton, Fourier 
transforms of p(x|0) and p(x, y|0) were expressed as sums 
of N and 0(N 2 ) terms in [2J], which allows an exact cal- 
culation of w(X). 

Useful limiting cases were also studied in [24[. For a 
broad external potential, with U(x) varying slowly over 
the width of the soliton, that can be estimated by the 
width £ of the classical soliton given in Eq. (f3T>)) . w(X) is 
essentially proportional to the square of the second order 
derivative of U, 



w{X)~C(N)^[U"(X)} 2 



(61) 



The coefficient C(N) only depends on N, it is given as a 
sum of 0(N 3 ) terms, and it has the large- N asymptotic 
behaviour 



C(N) 



N- 



N 



2vr 2 

— +4C(3) 



(62) 



where £ is the Riemann Zeta function. For a narrow 
external potential, centered at the origin of coordinates 
with a width much smaller than the soliton width £, one 
has for any N: 



w(X) ~ p(X\0) 



+ OO 



dxU 2 (x). 



(63) 



Finally, irrespective of the width of U(x), one can sim- 
plify Eq. (1601) in the large N limit by using an asymptotic 
expression for the pair distribution function |24j |: 



+00 



w(X) - 2NC dx dy U" (X + 2c£) U"(X + y£) 



2 + y - x 



(64) 



(ev + l)(e- x + 1) 
where Ng (and thus £) are kept fixed while N — > +00. 



C. General results for the error potential W(X) 

We now review some exact results derived in [24j on the 
function w(X) appearing in the numerator of the error 
potential W(X), see Eqs. (|47|48p . The function w(X) 
can be expressed in terms of the static structure factor 
of the ground state soliton for a fixed position X of its 
centre of mass: 

S(x,y\X) = ( P (x)p(y)) x (58) 
= S(x-y) P (x\X)+p(x,y\X) (59) 

where p(x,y\X) is the pair distribution function of the 
soliton of centre of mass localized in X. One has indeed 

w(X)= I dxdyU{x + X)U{y + X) 

x[S(x,y\0)-p(x\0)p(y\0)}. (60) 



D. An improved bracketing applicable to a Dirac 
external potential 

A limitation of the transmission and reflection coef- 
ficient bracketing of subsection 1111 Bl is that it becomes 
useless when the external potential U(x) is too narrow. 
For example, in the limiting case of a repulsive Dirac 
potential, U(x) = vS(x), v > 0, it is apparent that the 
quantity (V 2 )x in Eq. (|47p in infinite, since it contains 
a sum over all the particles of 8 2 (xi). As a consequence, 
the quantity e defined in Eq. (|51[) is +00 and the theorem 
applicability condition e < 1/2 is not satisfied. 

Here, we show that a slight improvement of the deriva- 
tion allows to remove this limitation. The resulting 
bracketing is thus more stringent, the price to pay be- 
ing that the new upper bound on \t — to\ is more difficult 
to evaluate in practice. 



7 



One simply uses the fact that 
QHQ-EQ>QH iut Q + NmfU(x)Q-EQ (65) 

X 

where as usual, for two Hcrmitian operators A and B, 
A > B means that the operator A — B is non-negative, 
that is (u\(A - B)\u) > for any ket \u). Eq. jfej) re- 
sults from the fact that the centre-of-mass kinetic energy 
operator P 2 /2M is non-negative, and that each oper- 
ator U(xi) is larger than or equal to m.i x U(x). Since 
we still impose Eq. (|22|) on the energy E, and since 
QHimQ > [E {N) + A]Q, the operator in the right- 
hand side of Eq. ()65[) is positive (within the subspace 
over which Q projects). 

For two positive Hermitian operators A and B such 
that A > B, one has that [35[ 



B~" > A- 



Applying this relation for A and B being the left-hand 
side and the right-hand side operators in Eq. (|55|). con- 
sidered within the subspace over which Q projects, one 
finds that 



o < g < 



Q 



NM x U(x)Q+QH int Q-EQ 



(67) 



where Q is defined in Eq. (|43[) . From the Cauchy-Schwarz 
inequality, writing Q = (G 1 ^ 2 ) , we have for arbitrary 
kets 1^1,2) such that Q 1 ^ 2 \ui 2 ) are normalisable: 



\{ui\g\u 2 )\ 2 < («i|e|ui>(u 2 |ei«2>. 



(68) 



We apply this inequality to the kets \u\) and {112) defined 
in Eqs. (|40|41[) . that do not need to be normalisable, and 
we use the upper bound on the operator Q to obtain 
for arbitrary centre-of-mass wavefunctions <E>i 2 (X) (not 
diverging too fast at infinity): 



where the improved error potential is positive: 
W imp (X) = (cj>\VQ 



Nmi x U{x)Q + QH int Q - EQ 



1/2 

5 

(69) 

QV\<t>). 

Thanks to the occurrence of the internal kinetic energy 
operator of the particles within H{ nt in the denominator, 
this error potential remains finite even when the barrier 
is a repulsive S potential. 

The reasoning of subsection IIII Al may then be repro- 
duced, replacing the error potential W by the improved 
one. Similarly to Eqs. (|51l52p we thus define 



c imp 



^?imp 



h 2 K\t \ 
(j»o|W imp (X)|l» ) 
(*ol^ imp (X)|$ ) 



(71) 
(72) 



where, as in subsection IIII Al <l>o(^0 and §>o(X) are the 
centre-of-mass scattering wavefunctions with incoming 
wavevector K and —K respectively, for the potential 
V{X). One then has: 

1/2 

Improved theorem: ifeimp^imp < 1/2 then 



1*1 - Ml < I* — *o| < 



\r\ - \r \\ < \r-r \ < 



U I 1/2 
iHllEimp^imp 

1 _ 2e n 1/2 

Kokimp 
Z£ imp']imp 



1 



(73) 

(74) 



where r, t are the exact reflection and transmission co- 
efficients, and ro,to are the reflection and transmission 
coefficients for $0, that is for the potential V(X). As in 
subsection IIII Al a simpler form is obtained for an even 
(66) external potential U (x) = U(—x), in which case ?7i mp = 1. 



E. General results for the improved error potential 

VFimp(X) 

A general calculation of Wi mp (X) with the Bethe 
ansatz, amounting to evaluating an internal dynamic 
structure factor of the ground state soliton with fixed 
centre-of-mass position, may be doable with the tech- 
niques developed in [H, [l6| but this is beyond the scope 
of this paper. On the contrary, a large N limit (for fixed 
Ng and £) is straightforward to obtain from the Bogoli- 
ubov technique exposed in Sec.[TV] In Eq. (|90|) one simply 
has to omit the centre-of-mass kinetic energy term and 
to replace V(X) by the lower bound N mf x U(x): 



W imp (X) 



N- 



dk 



[r*(AT)| s 



+™ U 2^ iVinf x U{x) + 4^ + A 



h 2 K 2 

2M 

(75) 
and A ~ 



where the amplitude Tk(X) is given by 
h 2 /(8m£ 2 ), see ®. 

When the external potential is a repulsive delta po- 
tential, U(x) = vS(x), v > 0, with mf x U(x) = 0, the 
integral over k can be calculated: 



Wi mp (X) 



2Nmv 2 £ 
-fh 2 



<t>l(X) 



1 32(7 + 2) ^^^; 

1 - -^T^C 00 (X) 



(7 + l) s 



where 7 £ (0, 1) is such that 
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1 



h 2 K 2 

2MA' 



(76) 
(77) 



As expected, (|76|) diverges when the incoming centre-of- 
mass kinetic energy tends to the gap A, but it diverges 
as I/7, whereas the error potential W(X) generically di- 
verges as I/7 2 for a non- negative U(x), see Eq. (|48|) . 

This indicates that the improved bound can have some 
interest also for a broad barrier. In the large N limit, 
when U(X) has a width b ^> £, we find using (I8T)1) and 
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assuming for simplicity that inf^ U(x) = 0: 



W imp (X) 



&>C A 
dk, 



■{U"{X)f 



x / ^[(l + A; 2 ) 2 ( 7 2 + fc 2 )cosh 2 ( 7 rfc/2)]- 1 , (78) 
Jr 27r 

where the integral may be expressed analytically if nec- 
essary, in particular in terms of the derivative of the 
digamma function, using the residue theorem. 



IV. LARGE TV LIMIT OF THE EFFECTIVE 
POTENTIAL FOR A 0(1/N) BARRIER 

We calculate a large N expansion of the non-local part 
SV of the effective potential in Eq. (f2l))) . using Bogoliubov 
theory, in the case where the external potential U expe- 
rienced by each particle scales as 1/N and the soliton 
width £ is fixed (because Ng is fixed). This physically 
convenient scaling with N ensures that the potential V 
has a well-defined non-zero limit for N — > +oo. 



A. Bogoliubov theory in brief 

We use Bogoliubov theory to dress with quantum fluc- 
tuations the classical soliton of single particle wavefunc- 
tion (|38|) . Since <fio is centered at the origin of coor- 
dinates, we shall shift the positions of the particles as 
Xi — > Xi + X where X is the fixed position of the centre 
of mass of the quantum soliton. In the number conserv- 
ing theory [2^, [26[ , one splits the bosonic field operator 
as 



i>{x) = a Q (j> (x) + 4>x(x) 



(79) 



where do annihilates a particle in the mode 4>o an d the 
field i(jj_(x) is orthogonal to the field cj>p(x). We intro- 
duce the modulus-phase representation [23, [28| , which is 
an excellent approximation for large N (when the prob- 
ability of having an empty mode </>o is negligible) : 



a 



(80) 



where the Hermitian phase operator 6 is conjugate to the 
number operator ho — ajao, [ri.o,0] = *■ The phase 9 is 
formally eliminated by its inclusion with the field ip± in 
the number conserving field 



k(x) = e- l0 4> ± {x). 



(81) 



Conservation of the total number of particles allows to 
eliminate no in terms of A and of the total number op- 
erator N. In the large N limit (with Ng fixed), it is 
found that A = (^(l/t; 1 / 2 ) whereas do scales as N 1 / 2 , 
which allows a systematic expansion of the Hamiltonian 



in powers of A. Keeping terms up to order <3(A 2 ) in Ho 
leads to the Bogoliubov approximation for N particles, 



P 2 



2Nm 



R 



dk 

2^ 



e k b{b k (82) 



with the gapped Bogoliubov spectrum in terms of the 
quasiparticle wavevector k and the Gross-Pitaevskii 
chemical potential /iq- 



I Mo I 



h 2 k 2 

2m 



where uq 



(83) 



The quasi-particle annihilation and creation operators bk 
and bl obey the usual bosonic commutation relations on 
the open line, [Sfcj&j,,] = 2n5(k — k'). Due to the trans- 
lational symmetry breaking, a Goldstone mode appears, 
with a massive term tx P 2 in the Hamiltonian, the field 
variable P (scaling as TV 1 / 2 ) representing at the Bogoli- 
ubov level the total momentum of the system, and be- 
ing conjugate to the field variable Q (scaling as 1/N 1 / 2 ) 
giving at the Bogoliubov level the fluctuations of the 
centre-of-mass position of the system: This reproduces 
the structure of Eq. ([4]). The modal field expansion is 
then 



k(x) = -N^^Q + j-^xM^P 



R 



dk 

2^ 



u k (x)b k + v* k (x)b\ 



(84) 



The Bogoliubov mode functions are known exactly [29( 
and are given with the present notations in [24j|. They 
are orthogonal to ipo, and one has also that u k + v k is 
orthogonal to x(f>o(x). This is apparent on the useful 
form: 



(/) {x)[u k (x) + v k (x)} 



4£ 2 



{l + 2i\k\£) 2 dx 2 



[e lkx Ux)\ ■ 
(85) 



B. Bogoliubov expression of SV 

To calculate SV given by Eq. (|3"T|) . we first have to 
express the operator V in the Bogoliubov framework. To 
the same level of approximation as for the Hamiltonian 
Hq, that is neglecting terms that are cubic or more in A, 
we obtain 

y2u(X+ Xi ) ~ (N- [ dxk^A) [ dxcj)l(x)U(X + x) 
i=1 Jr Jr 

+ N 1/2 ( dx[k{x)+k\x)](j)o{x)U{X + x) 
Jr 

+ ( dxU{X + x)k\x)k{x). (86) 



It is clear that, in the Bogoliubov framework, applying 
the projector V of the full quantum theory amounts to 
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projecting onto the vacuum |0)sog of all the quasiparticle 
annihilation operators b k , so that applying Q amounts to 
projecting onto the subspace with at least one quasipar- 
ticle excitation. In the Bogoliubov evaluation of QV\cf>), 
the leading term explicitly scaling as N in Eq. (1551) gives 
a vanishing contribution, so we keep the subleading N 1 ^ 2 
term of (|86|) . Similarly, due to the projector Q, we keep 
only the contributions involving b' k in Eq. (p?4"|) to obtain 



QV\4>) 



dk 
2^ 



T k (X)bl\0) Bog 



with the amplitudes depending parametrically on X: 



T k {X)=N 1 ' 2 / dx[u* k {x) +v* k (x)}4> (x)U{X + x) 

4^2^1/2 



(l-2t\k\0 2 



dxe- lkx Mx)U"(X + x) 



where we used (|85p . The integral over x is typically cut to 
|x| < £ by the rapidly decreasing function 4>o(x). For an 
external potential narrower than £, |rfc(X)| varies with 
X at the length scale £. For a broad external poten- 
tial, varying at a scale b 3> £, we expect that T k (X) 
varies with X at the length scale b. This can be made 
quantitative by expanding U"(X + x) in the integrand of 
Eq. (|88p to zeroth order in x: 



T k (X) 



47rf /2 iV 1/2 ;7"(X) 
6»C (1 - 2i|fc|0 2 cosh(7rfc^)' 



(89) 



It remains to estimate the denominator EQ — QHQ 
with Bogoliubov theory. Within the subspace with one 
Bogoliubov excitation, the leading term explicitly scaling 
as N in Eq. (|86j) gives a non-zero contribution which is 
actually O(l) since U = 0(1/N). This contribution does 
not affect the quasiparticle, it is a scalar depending on X 
only, and it simply corresponds to the mean-field (Gross- 
Pitaevskii) approximation for V(X). The subleading 
term in Eq. (|86[) . which changes the number of quasipar- 
ticles by ±1, is more involved since it couples the single- 
quasiparticle subspace to the two-quasiparticle subspace. 
However, its contribution is OiN 1 / 2 !/) = Oil/N 1 / 2 ) and 
may be neglected at the present order. 

We finally obtain in the large N limit, for U = Oil/N) 
and fixed Ng, the leading term of SV, scaling as 1/N: 



SV 



dk 
2^ 



— +E (N) + e k + V(X) 



r k (x). 

(90) 



Here P = —ihdx is the centre-of-mass momentum oper- 
ator of the full quantum theory. The expansion of V up 
to the same order (neglecting o(l/N)) is directly given 
by (|39p . One may then solve numerically the resulting 
approximate form of Eq. (1261) . which is made delicate 



by the non-local nature of (1901 . More simply, one may 
treat (|9"0"f as a first order perturbation in the V scatter- 
ing problem using the formulation of the Appendix lAl [see 
below Eq. (|A6[) ]. to obtain: 



t-ta 
r-r 



w <*SI*v|*<,) 

iM 

Wk 



:*5i*vi*o>, 



(91) 
(92) 



the ket |$q) corresponding to the wavefunction $q(X). 



V. BORN-OPPENHEIMER-LIKE APPROACH 

In molecular physics, one often uses the so-called Born- 
Oppenheimer approximation: One diagonalises the elec- 
tronic problem for fixed positions of the nuclei, obtaining 
a ground state electronic energy that then serves as a po- 
tential for the nuclei [l9|, [3(| • It is natural to try to apply 
a similar approach to our soliton scattering problem. The 
"heavy" particle then corresponds to the centre of mass 
of the soliton, and the "light" particles are the internal 
degrees of freedom of the soliton. We thus split the N- 
body Hamiltonian as H = iJ hcav y + H lisht , with 



^hcavy = +NU(X) 

2M v ' 



N 



H^ ht = H int + Y}U{x t )-U{X)]. 



(93) 



(94) 



In H heavy , we have included the external potential that 
the soliton would feel if all the particles were localized 
in the centre-of-mass position X. We expect this term 
NU(X) to constitute already a good approximation when 
the external potential is much broader than the soliton 
width £. To go beyond this zeroth order approximation, 
the idea is to calculate the ground state energy i5g lght (X) 
of H h&ht for a fixed value of the centre-of-mass position. 
This energy then provides a correction to the potential 
NU(X) experienced by the centre of mass of the soliton. 
An important condition for the Born- Oppenheimer ap- 
proximation to hold is that i?Q lght (A) is well separated 
from the excited state energies of H hght for a fixed X , so 
that the presence of the solitonic internal gap again plays 
an important role here. 

More precisely, we call |x(A)) the ground state of 
the internal Hamiltonian H llsht corresponding to the 
eigenvalue i?g lght (AT), and we put forward the Born- 
Oppcnheimer-like ansatz for the A^-body state vector: 



(X|V)bo = HX)\x(x)) 



(95) 



where the ket \X) represents the centre of mass perfectly 
localized in X and the internal ket \y(X)), normalized to 
unity, parametrically depends on X [361 ] . We then insert 
the ansatz into Schrodinger's equation E\ip) = H\ip) and 
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project with (\(X)\ to obtain [37 



2M dX 2 

h 2 



NU(X) 
d 2 



E^ ht (X) 



T7T7(X(X)\--\X(X)} $(X). (96) 



2M 



'dX 



Note that we keep here the so-called Born-Oppenhcimer 
diagonal correction coming from the X dependence of the 
internal state in the ansatz. 

In practice, to evaluate EQ Sht (X) and the correspond- 
ing eigenvector, we use perturbation theory, treating 
J2iLi[U( x i) — U(X)} as a perturbation of H m f For ex- 
ample, to second order in this perturbation, we obtain: 



E^{X) ~ Eo(N) - NU(X) + <0|V|0) 



\VQ 



Q 



E (N)Q-QH int Q 



QV\4), (97) 



where, as in the previous sections, the internal ket \<j>) 
is the free space ground state soliton of energy Eo(N), 
V projects orthogonally onto \<f>) and the supplementary 
projector Q = 1 — V projects onto the internal excited 
states of the system. Remarkably, the third term in the 
right-hand side of Eq. (l97l) exactly coincides with V(X), 
To first order in the perturbation theory, the internal 
ground state in presence of the external potential is 



\ X {X))^T{X) 



Q 



E Q (N)Q-QH int Q 



QV\<f>) 



(98) 



The normalization factor J- should for consistency only 
weakly deviate from unity, which imposes a limit on the 
strength of the external potential. 

We apply the above approximation scheme in the large 
TV limit, where it makes the most sense, fixing Ng (and 
thus the soliton width £). We can then use the Bogoli- 
ubov approach, and following the lines of section IIVI 



E, 



light 
o 



(X) ~ B Bog (JV) - NU(X) + V(X) 

dk \T k (X)\ 2 



2tt 



(99) 



where Tk{X) is given by Eq. (|85]l. Also, the Born- 
Oppenheimer diagonal correction is approximated with 
Bogoliubov theory as 



2M {x{X) ^dX 2 



X(X)) 



2M 



dk \^r k (x)\ 2 



2tt 



(100) 



It is about N times smaller than the Bogoliubov term 
appearing in Eq S t (X) ) see the last term in Eq. and 
we neglect it in the large N limit. To summarize, in the 
Born-Oppcnhcimcr approximation for large N, we find 
for the soliton wavefunction an equation of the form (|26p 
with the non-local SV approximated by the local form 
[after use of Eq. 



<5V BO 



dk |r fc (X)| ; 



2tt 



-Ck 



8A 



dxdy Mx)Mv)U"(X + x)U"{X + y)e-W*[(x - y) 2 + 6^\x -y\ + 12£ 2 ] 

(101) 

I 



The integral can be evaluated for a delta external poten- 
tial, U(x) — vS(x), leading to (38J: 



SV 



BO 



(102) 



It is interesting to compare Eq. (|101j) to the result 
Eq. (|90l) that was obtained in a different way. At first 
sight, Eq. (jnOI) and Eq. (fTUTj) look widely different, be- 
cause of the more complicated energy denominator in 
Eq. (JOT]) that involves both P 2 /2M and V(X), rather 
then a simple c- number quantity such as We have 
identified a limiting case where the two expressions are 
close, when the typical wavevector K of 3>(X) is much 
larger than l/£, and SV has a small perturbative effect on 
the scattering state $. Since the relevant k are 0(l/£) in 
the integral over k, Tk(X) varies at a length scale of order 



£ or larger, see discussion below Eq. (|88p . whereas <&(X) 
varies over a much smaller length scale 1/K. The spatial 
derivatives of Tk(X)<$>(X) are thus well approximated by 
taking the derivatives of $>(X) only, e.g. 



PT k (X)\$) 



r fe (x)p|$>, 



(103) 



where P = —ihdx is the centre-of-mass momentum op- 
erator. We can then approximately commute the en- 
ergy denominator with Tk(X) in Eq. The last step 
is to realise that, if £V is small enough, $>(X) will be 
close to the scattering state of energy E for the potential 
Eq(N) + V(X) so that 



E — 



^L+Eo(N)+e k + V(X) 



1$) 



~l*> 
(104) 
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and we recover the Born-Oppenheimer result Eq. (|101l) . 



VI. APPLICATIONS OF THE FORMALISM 

A. Centre-of-mass wavepacket splitting 

We apply our formalism to the proposal of [6j for the 
production of Schrodinger's cat-like states by elastic scat- 
tering of a soliton on a barrier: One sends the ground 
state soliton with a quasi-monochromatic centre-of-mass 
wavepacket, that is centered in A'-space around K with 
a width AK <C K, on a barrier of adjusted height such 
that the transmission and reflection amplitudes have the 
same modulus 1 / y/2. This prepares the gas in a coherent 
superposition of all the particles being to the right and 
to the left of the barrier with equal probability ampli- 
tudes. For the experimental decoherence rate estimated 
in this in principle allows to prepare a gas of N ~ 100 
lithium 7 atoms in a coherent superposition of being at 
two different locations separated by ~ 100 /im. We thus 
restrict to the most interesting large- N limit. As the po- 
tential barrier may be produced with a focused Gaussian 
laser beam, we can assume that U(x) is a repulsive Gaus- 
sian of width b (the so-called waist of the laser beam) : 



U(x) = U e- 2x2 ' b \ U >0, 



(105) 



so that the potential V(X) is also even and bell shaped. 
We shall also assume, as in @, that 



h 2 K 2 _ A 
2M ~ ~2' 



(106) 



so that in the large N limit, 



K 



jVV2 



(107) 



Having a significantly smaller K would indeed uselessly 
slow down the Schrodinger's cat formation process. Hav- 
ing a significantly larger K is forbidden by the elasticity 
condition (l2"2"j) . 

Case t » ^: This broad barrier case is experimen- 
tally the typical one, since the waist of a focused laser 
beam is a few microns, whereas £ < |6j. Then 

V{X) ~ NU(X), and the width of V is also of order 
b. Since Kb oc N 1 ^ 2 b/^ is much larger than unity, the 
scattering problem of the centre-of-mass on V(X) is in 
the semiclassical regime, where the transmission and re- 
flection amplitudes at incoming wavevector K have the 



approximate expressions, see Eqs. (3.49,3.58,4.23) in [3l| : 

e -\(K) e i[<x(K)-0(K)] 



t 



ro 



[1 



-% e 



-2A(if)ll/2 
i[a(K)-p(K)] 



X(K) 

a(K) 
f3(K) 

K(X) 



[1 + e 



(108) 
(109) 
(110) 



2A(if)jl/2 

dX (-i)K(X) > 
2KX- + 2 f dX[K{X)-K] (111) 

J — oo 



X(K) 



In 



X(K) 



2 2MV(X) 



-,1/2 



(113) 



Here X_ < and X + = —X^ are the two classical turn- 
ing points for an incoming energy below the maximum 
Vo of V, and K (X) is either in R+ or in iR+ if X is in the 
classically allowed or forbidden region. For an incoming 
energy larger than Vq, one has to use analytic continua- 
tion [3l| . The phase P(K) is given in 31 1, and the extra 
phase a(K) is due to our different choice for the phase 
reference point. We also recall the WKB approximation 
for &q(X) in the classically allowed region: 



Mx) x< ~-[k/k{x)Y' 2 



e iKX e if^ dx[K(x)-K\ 



. rQe -iKX e -iJ^dx[K(x)-K] 



(114) 



^ (X) X> ^ + [K/K(X)}^ 2 t Q e lKX e i ^ dx[K{x) ~ K \ll5) 

We then see that a transmission probability of 1/2 is 
achieved in the semiclassical formula (|108|) for an energy 
equal to Vo, that is for a momentum 



/2 



(2MV ) 1/2 



(116) 



Away from this value of K, \to\ 2 will drop rapidly to 
zero or rise rapidly to one. A local formula is obtained 
by approximating the top of V(X) around X = by a 
parabola, so that 



M 2 



i 



exp[(X 1/2 - K)/8K] 



with 



SK 



1 

2^ 



\V"( 



2V(0) 



1/2 



(117) 



(118) 



For large N, one finds SK ~ 1/(V2tt&) so that one is ex- 
perimentally in the regime AK 3> SK . In what follows, 
we thus adjust the barrier height to have K\i% = K. 
Then, with Eq. (flOo) . 



NU ~ Vo 



A 
~2' 



(119) 
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Finally, we evaluate the parameter e of (|5 1 1) appearing 
in the bracketing (|56I57I) (that can be used here since 
U(x] and V(X) are even), for the physically most rele- 
vant case K = K = Ki/2- We use the large- TV estimate 
for a broad barrier, see Eqs. (161l62p . and the simplest 
estimate (neglecting rapidly oscillating terms): 



(120) 



with K(X) defined in (| 1 13[) and a proportionality factor 
equal to 1 + |r | 2 = 3/2 for X < and to |t | 2 = 1/2 
for X > 0. Since K(X) vanishes linearly in the classical 
turning point X = 0, we get a logarithmic divergence 
in the resulting approximation for ($o|W|^o), that we 
cut by introducing the quantum length scale ctho -C b 
associated to the Schrodinger's equation in the inverted 
parabola approximating V(X) close to its maximum [39| : 



h 1 ' 2 



( M| V" (0)|) 



_ ( _^_\ 1/2 ( 2bg \ 1/2 
1 ' i ~ \V2KJ ~ \NV 2 ) 



. (121) 

Keeping only the logarithmically diverging contribution 
amounts to approximating the matrix element as 

($ \W(X)\<f> )~2W(0)^ dX^-j=, (122) 

which leads to the estimate 



2^ 3 (2^ 1/2 

N 



¥ 



2C(3) 



ln(6/a ho ). (123) 



In conclusion, for |to| 2 — 1/2 in the broad barrier case, 
we find in the large N limit (where e< 1): 



\t-t \=O 



m 3 

iVi/a 



In 



(124) 



as already given in [6|. 

Case b <C For a narrow barrier, 



V(X) ~ p(A|0) 



(125) 



In the large N limit, replacing p(X\0) by its classical field 
approximation, that is the leading term in the right-hand 
side of (j3"9")l . gives 



V(X) 



Vn 



with 



cosr/(X/& e ff) 



6 cff - 2£. 



(126) 

(127) 
(128) 



Although the resulting scattering problem for $> (X) 
then becomes exactly solvable [32[, we simply reuse the 



semiclassical reasoning of the previous (broad-barrier) 
case, since Kb e s ~ (N/2) 1 / 2 ^ 1 is again in the semi- 
classical regime. At half transmission probability, 



NUp 
A 



1/2 



(129) 



so that NUo is now much larger than the gap A, con- 
trarily to the broad barrier case. The harmonic oscilla- 
tor length used in the cutting procedure is found to be 
a ho ^ b cS (2/Ny/ 4 . We use the equivalent of Eq. ([122]) 
and we estimate W(0) from Eq. (|63|) . At half transmis- 
sion probability, we finally obtain from the simple brack- 
eting (|56I57I) : 



\t-t \ < 



Zjb HN/2) 

N l/2 ( 27r )l/2 ' 



(130) 



As expected, this bound diverges for b — > (at fixed 
N), since U(x) then approaches a Dirac potential, for 
which the use of the improved bracketing (|73l74p is more 
appropriate and leads at half transmission for large N to 



I* — *o| 



< 



ln(7V/2) 
8iVV2 ' 



(131) 



where Eq. (|76|) was used with v = Uob (it/2) 1 / 2 . 
The Born-Oppenheimer prediction: For the delta exter- 
nal potential U(x) — v8(x), it is interesting to compare 
the upper bound (|131l) to the result of Sec. [V] In the 
present large N limit, one can treat Vbo with first order 
perturbation theory similarly to Eqs. (|91l92p and one can 
use the expression (|102l) for Vbo- In the resulting per- 
turbative expression for t — to, one can approximate the 
various quantities by their N — > +oo limit, in particular 
the scattering wavefunction <&o(X) at incoming wavevec- 
tor K may be replaced by the scattering wavefunction 
<$> { °\X) of the 1/cosh 2 potential of Eq. (fl"2"6"l) . which is 
exactly expressed in terms of an hypergeometric function 
(32| . with the transmission amplitude 



M _ 
L o — 



ra + is-iKb ea )rn 



iKb r 



T(l - iKb cS )T(~iKb cS ) 



with 



( 2Mb 2 cg 1 



1/2 



(132) 



(133) 



Here, to zeroth order in 1/iV, we have Kb s g ~ (N/2) 1 / 2 
as in Eq. (|107p . and s ~ Kb c g due to the half- 
transmission condition | io°^ I — 1*0 1 = l/v^, so that 
Nmv 2 /h 2 ~ h 2 /(16M£ 2 ). Further expressing Eq. (fT02|) 
in terms of the variable 9 = tanh(A/6 c ff), we obtain 



t-to 



8(JV/2)i/2 



de(i-"-e 2 +^)^\x)^>(-x). 

(134) 
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The terms proportional to 9 2 and 9 A vanish in 6 = 0, 
which allows to directly use the WKB forms (| 1 1411 15p : 
at half-transmission, K(X) = K\6\ so that 



*$\x)&»>{-X) 



(0), 



+(°) / 

l [ i l_ e -iJf6««ln(l-6 

y/2 



(135) 

and one finds as expected that the rapidly oscillation bit 
gives a negligible contribution to the integral. Note that 
the semiclassical approximation gives 



(0) 



-2iKb e a In 2 



V2 



(136) 



in agreement with the Stirling asymptotic equivalent of 
(fT32)) for s = Kb cS -> +oo. 

On the contrary, the constant term in between the 
parenthesis in Eq. (|134p cannot be treated with the sim- 
ple WKB approximation (|135|) : As already discussed 
above, this simple approximation is inaccurate over the 
interval \X\ < a no , where it would incorrectly lead to 
a logarithmic divergent integral. As a straightforward 
alternative to more elaborate semiclassical methods, we 
can use the fact here that an infinitesimal change 5V$ of 
the amplitude Vq of the 1 / cosh 2 potential will lead to 

change of that may either be evaluated by perturba- 
tion theory, or by taking the derivative of Eq. (|132[) with 
respect to Vq, that is with respect to s. This leads to the 
exact relation 

f de^\x)^\-x) = -^-ltf[^+is-iKb ei! ) 

J -l s I 

- - is - iKb eS )\ (137) 

where the digamma function i/)(z) behaves as lnz + o(l) 
for \z\ — > +oo. 

In conclusion, for the soliton scattering at a centre-of- 
mass kinetic energy A/2 on a delta external potential 
such | to | = l/\/2, the Born-Oppcnhcimcr-like approach 
predicts that, in the large N limit with Ng fixed, 



t-tn 



t 8(JV/2)Va 



ln(2JV) - V(l/2) -i\ 



3 

2 4 

(138) 

where we recall that -0(1/2) = —2 In 2 — C and C = 
0.57721 ... is Euler's constant. This is compatible with 
the bound (|13ip . Since the equivalence conditions of 
the Born-Oppenheimer-like approach with the system- 
atic Bogoliubov approach of Sec. lIVI are here satisfied, as 
discussed in the paragraph below Eq. (|102j) . the typical 
centre-of-mass wavevector diverging as A 1 / 2 , the result 
f|138[) is asymptotically exact |40] |. 



is for the scattering of a dimer, on a delta barrier U(x) = 
v5(x), v > 0. In this case, the set of internal coordinates 
reduce to the relative coordinate x — — X\ of the two 
particles, with x\ = X — x/2 and X2 = X + x/2. The 
internal Hamiltonian is simply 

H ini = -^ +9 S(x). (139) 

Its normalized ground state wavefunction is <fi(x) = 
Q^ 2 e~ go I 35 1 , with an energy Eq(2) = —h 2 q 2 /m in agree- 



% 

ment with Eq. ([S]), where we have set 



Qo = 



mg 
2ft 2 ' 



(140) 



This immediately gives the mean potential 

V(X) = 4v\(j)(2X)\ 2 = 4vq e- 4q °W. (141) 

Since the continuous spectrum of H lnt starts at zero en- 
ergy, one has the gap A = —E (2), in agreement with 
Eq. ([S]). The Green's function of H- lnt at energy 



fl 2 q 2 

E = , 0<q<q , 



m 



(142) 



is also easily calculated from the differential equation 

[E + -dl - gS(x)](x\(E - H^ly) = 6(x - y) (143) 
m 

with the boundary conditions that it does not diverge 
exponentially for \x\ —> +oo. E.g. for y < 0, one simply 
has to integrate the differential equation over x over the 
intervals x < y, y < x < 0, < x, where the general 
solution is the sum of two exponential functions of x, 
and then match the solutions in x = y and x = 0, using 
the continuity of the Green's function with x, and the 
discontinuity of its first order derivative with respect to 
x as imposed by the Dirac terms. To finally obtain the 
matrix elements of Q/(EQ — QH lrA Q) in position space, 
one simply has to remove from the Green's function of 
Hint the contribution cf)(x)(f)(y)/[E — Eq(2)] of the ground 
state of Hi nt . We finally obtain 



w-q l qo - q 

Jqqo_ igo \x\ 

~ 2 2 

<7o - q 



(144) 



A numerical or an analytical solution [4l[ of the scat- 
tering problem for <f> (X) can then be combined to this 
expression for Wi mp (X), to obtain explicit numbers for 
the improved bracketing (|73I74[) , 



B. Application of the improved bracketing to N = 2 



VII. CONCLUSION 



Explicit calculations of the improved error potential 
Wi m p (A) of Eq. ([701) ma Y be performed for N = 2, that 



We have considered the scattering of a one-dimensional 
quantum soliton (the bound state of N attractive-^ 
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bosons) on a potential barrier in the elastic regime, 
where energy conservation prevents observation of soli- 
ton fragments at infinity. The scattering at a given in- 
coming centre-of-mass wavevector K is then character- 
ized by reflection and transmission amplitudes r and t 
for the soliton centre-of-mass wavefunction 3>(X), with 
|i| 2 + |r| 2 = l. 

In the simplest approximation, one assumes that &(X) 
simply sees a local potential V(X) obtained by averag- 
ing the single particle external potential U (x) over the 
particle density profile of the quantum soliton, leading 
to approximations tq and to for the amplitudes r and t. 
Rigorous upper bounds on the resulting errors \r— 7*0 1 and 
\t — to I are derived and are expressed in an operational 
form (distinguishing various limits of broad or narrow 
barrier, for any N or for large N). 

In an exact treatment, also giving a precise meaning 
to &(X), it is shown that an additional, non-local po- 
tential SV appears in an effective Schrodinger's equation 
for $. The large N leading behaviour for SV is obtained 
using Bogoliubov theory, and it is compared to a Born- 
Oppenheimer-like approach that treats the centre of mass 
of the system as the heavy particle. 

Finally, simple applications of the formalism are given, 
mainly in the context of the Schrodinger cat state pro- 
duction scheme considered in @, 0] ■ 
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Appendix A: Bounds on the error on t and 



In this Appendix we shall prove Eqs. (l54l55p . We 
rewrite Eq. (p?S|) in position space, 



h 2 K 2 
2M 



MX) 



2M 



$"(A) + V(X)$(X) + S(X) (Al) 



with the contribution that we shall treat formally as a 
source term, 



S(X) = (X\6V\&) 



(A2) 



Here &(X) obeys the boundary conditions (|14I15I) . and 
the fact that SV is hermitian leads to Irl 2 + \t\ 2 = 1 as 



expected [42|. To integrate formally Eq. (|A1|) we need 
two independent solutions of the corresponding homoge- 
neous equation. One is $q(X), i.e. the scattering solution 
for an incoming wave from the left (i.e. with a centre-of- 
mass wavevector K > 0). The other solution is conve- 
niently taken as the scattering solution for a wave incom- 
ing from the right (i.e. with a centre-of-mass wavevector 



(-K) < 0), denoted as $o(X). If the potential U(x) is 
even, we simply have &o(X) = <E>o( — X). In the general 
case, we take Eq. ([53| . Then one may check that 



* {X) ~ e~ lKX - r ^KX for x ^ +QO (A3) 

*0 



0>o(X) ~ t e- tKX for X -> -c 



(A4) 



Then, after formal integration with the method of vari- 
ation of constants and calculation of the Wronskian of 
$ (X) and $ (X), 

W(X) = <f>o(X)& (X)-& (X)$ (X) = -2iKt , (A5) 
we obtain 

<P(X) = [1 + A(X)]$ {X) + B{X)® (X) (A6) 

with 



A(X) = 



iM 
h 2 Kt J_ 



x 



dx<f> (x)S(x) (A7) 



B(X) = 



iM 
h 2 Ktc 



x 



dx$ (x)S(x) 



(A8) 



One may check that &(X) obeys the right boundary con- 
ditions dump with 



t = t [l + A(+oo)] 
r = r + t B(-oc). 



(A9) 
(A10) 



An upper bound of |j4(X)| is obtained from Eq. (|4T))) by 
taking 



$ 2 (x) = $(x), 



(All) 
(A12) 



where 9 is the Heaviside step function. Furthermore we 
use the fact that 



since W(X) is positive. Then 

\A{X)\<er,^\ 



(A13) 



(A14) 



where e is defined in Eq. (I5T1) . 77 is defined in Eq. ([52 
and 



($\W(X)\$) 



(A15) 



(*o|W(X)|$o) 
Similarly, taking <J>i(x) = 6{x — X)$q(x), one has 

\B(X)\ < ea 1/2 . (A16) 

The last step is to derive an upper bound for a. To this 
end, we derive an upper bound on |$(X)| from Eq. (|A6[) : 



|*P0| < l^oWlW^^^I^^I+lioWl]. (A17) 
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We square this relation, multiply by the positive quantity 
W(X), integrate over X, divide by ($ \W(X)\$ ) and 
use the Cauchy-Schwarz inequality 



+00 



dX\$o(X)\\$ {X)\W(X) <rfl 2 

/+00 
dX\$ Q (X)\ 2 W(X) (A18) 
-00 

to obtain 

a < l+4er/ 1/2 a 1/2 +4e 2 r/a = ( 1 + 2erj 1/2 a 1/2 ) 2 . (A19) 



Taking the square root leads to 

a 1 / 2 < l + 2e7 7 1 /2 a i/2_ 

For ery 1 / 2 < 1/2 we thus get 

1 



,1/2 < 



l-2e?7 1 /2- 



(A20) 



(A21) 



This inequality together with Eqs. (|A9IA10IA14IAlGl) . 
leads to Eqs. (|54l55[) . In the particular case of an even 
external potential U(x), W{X) is also even and one has 
r] = 1, which leads to the simpler relations Eqs. ()56|57|) . 
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